Antenna Season Report Notebook

Josh Dillon, Last Revised January 2022

This notebook examines an individual antenna's performance over a whole season. This notebook parses information from each nightly rtp_summarynotebook (as saved to .csvs) and builds a table describing antenna performance. It also reproduces per-antenna plots from each auto_metrics notebook pertinent to the specific antenna.

import os
from IPython.display import display, HTML
display(HTML("<style>.container { width:100% !important; }</style>"))

# If you want to run this notebook locally, copy the output of the next cell into the next line of this cell.
# antenna = "54"
# csv_folder = '/lustre/aoc/projects/hera/H5C/H5C_Notebooks/_rtp_summary_'
# auto_metrics_folder = '/lustre/aoc/projects/hera/H5C/H5C_Notebooks/auto_metrics_inspect'
# os.environ["ANTENNA"] = antenna
# os.environ["CSV_FOLDER"] = csv_folder
# os.environ["AUTO_METRICS_FOLDER"] = auto_metrics_folder

# Use environment variables to figure out path to the csvs and auto_metrics
antenna = os.environ["ANTENNA"]
csv_folder = os.environ["CSV_FOLDER"]
auto_metrics_folder = os.environ["AUTO_METRICS_FOLDER"]
print(f'antenna = "{antenna}"')
print(f'csv_folder = "{csv_folder}"')
print(f'auto_metrics_folder = "{auto_metrics_folder}"')

antenna = "137"
csv_folder = "/home/obs/src/H5C_Notebooks/_rtp_summary_"
auto_metrics_folder = "/home/obs/src/H5C_Notebooks/auto_metrics_inspect"

display(HTML(f'<h1 style=font-size:50px><u>Antenna {antenna} Report</u><p></p></h1>'))

Antenna 137 Report

import numpy as np
import pandas as pd
pd.set_option('display.max_rows', 1000)
import glob
import re
from hera_notebook_templates.utils import status_colors, Antenna

# load csvs and auto_metrics htmls in reverse chronological order
csvs = sorted(glob.glob(os.path.join(csv_folder, '*.csv')))[::-1]
print(f'Found {len(csvs)} csvs in {csv_folder}')
auto_metric_htmls = sorted(glob.glob(auto_metrics_folder + '/auto_metrics_inspect_*.html'))[::-1]
print(f'Found {len(auto_metric_htmls)} auto_metrics notebooks in {auto_metrics_folder}')

Found 237 csvs in /home/obs/src/H5C_Notebooks/_rtp_summary_
Found 224 auto_metrics notebooks in /home/obs/src/H5C_Notebooks/auto_metrics_inspect

# Per-season options
mean_round_modz_cut = 4
dead_cut = 0.4
crossed_cut = 0.0

def jd_to_summary_url(jd):
    return f'https://htmlpreview.github.io/?https://github.com/HERA-Team/H5C_Notebooks/blob/main/_rtp_summary_/rtp_summary_{jd}.html'

def jd_to_auto_metrics_url(jd):
    return f'https://htmlpreview.github.io/?https://github.com/HERA-Team/H5C_Notebooks/blob/main/auto_metrics_inspect/auto_metrics_inspect_{jd}.html'

Load relevant info from summary CSVs

this_antenna = None
jds = []

# parse information about antennas and nodes
for csv in csvs:
    df = pd.read_csv(csv)
    for n in range(len(df)):
        # Add this day to the antenna
        row = df.loc[n]
        antnum = row['Ant']
        if antnum != int(antenna):
            continue
            
        if this_antenna is None:
            this_antenna = Antenna(row['Ant'], row['Node'])
        jd = [int(s) for s in re.split('_|\.', csv) if s.isdigit()][-1]
        jds.append(jd)
        this_antenna.add_day(jd, row)
        break

---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
~/anaconda/envs/RTP/lib/python3.7/site-packages/pandas/core/indexes/base.py in get_loc(self, key, method, tolerance)
   3360             try:
-> 3361                 return self._engine.get_loc(casted_key)
   3362             except KeyError as err:

~/anaconda/envs/RTP/lib/python3.7/site-packages/pandas/_libs/index.pyx in pandas._libs.index.IndexEngine.get_loc()

~/anaconda/envs/RTP/lib/python3.7/site-packages/pandas/_libs/index.pyx in pandas._libs.index.IndexEngine.get_loc()

pandas/_libs/hashtable_class_helper.pxi in pandas._libs.hashtable.PyObjectHashTable.get_item()

pandas/_libs/hashtable_class_helper.pxi in pandas._libs.hashtable.PyObjectHashTable.get_item()

KeyError: 'Ant'

The above exception was the direct cause of the following exception:

KeyError                                  Traceback (most recent call last)
/tmp/ipykernel_4136/3863632283.py in <module>
      8         # Add this day to the antenna
      9         row = df.loc[n]
---> 10         antnum = row['Ant']
     11         if antnum != int(antenna):
     12             continue

~/anaconda/envs/RTP/lib/python3.7/site-packages/pandas/core/series.py in __getitem__(self, key)
    940 
    941         elif key_is_scalar:
--> 942             return self._get_value(key)
    943 
    944         if is_hashable(key):

~/anaconda/envs/RTP/lib/python3.7/site-packages/pandas/core/series.py in _get_value(self, label, takeable)
   1049 
   1050         # Similar to Index.get_value, but we do not fall back to positional
-> 1051         loc = self.index.get_loc(label)
   1052         return self.index._get_values_for_loc(self, loc, label)
   1053 

~/anaconda/envs/RTP/lib/python3.7/site-packages/pandas/core/indexes/base.py in get_loc(self, key, method, tolerance)
   3361                 return self._engine.get_loc(casted_key)
   3362             except KeyError as err:
-> 3363                 raise KeyError(key) from err
   3364 
   3365         if is_scalar(key) and isna(key) and not self.hasnans:

KeyError: 'Ant'

# build dataframe
to_show = {'JDs': [f'<a href="{jd_to_summary_url(jd)}" target="_blank">{jd}</a>' for jd in jds]}
to_show['A Priori Status'] = [this_antenna.statuses[jd] for jd in jds]

df = pd.DataFrame(to_show)

# create bar chart columns for flagging percentages:
bar_cols = {}
bar_cols['Auto Metrics Flags'] = [this_antenna.auto_flags[jd] for jd in jds]
bar_cols[f'Dead Fraction in Ant Metrics (Jee)'] = [this_antenna.dead_flags_Jee[jd] for jd in jds]
bar_cols[f'Dead Fraction in Ant Metrics (Jnn)'] = [this_antenna.dead_flags_Jnn[jd] for jd in jds]
bar_cols['Crossed Fraction in Ant Metrics'] = [this_antenna.crossed_flags[jd] for jd in jds]
bar_cols['Flag Fraction Before Redcal'] = [this_antenna.flags_before_redcal[jd] for jd in jds]
bar_cols['Flagged By Redcal chi^2 Fraction'] = [this_antenna.redcal_flags[jd] for jd in jds]
for col in bar_cols:
    df[col] = bar_cols[col]

z_score_cols = {}
z_score_cols['ee Shape Modified Z-Score'] = [this_antenna.ee_shape_zs[jd] for jd in jds]
z_score_cols['nn Shape Modified Z-Score'] = [this_antenna.nn_shape_zs[jd] for jd in jds]
z_score_cols['ee Power Modified Z-Score'] = [this_antenna.ee_power_zs[jd] for jd in jds]
z_score_cols['nn Power Modified Z-Score'] = [this_antenna.nn_power_zs[jd] for jd in jds]
z_score_cols['ee Temporal Variability Modified Z-Score'] = [this_antenna.ee_temp_var_zs[jd] for jd in jds]
z_score_cols['nn Temporal Variability Modified Z-Score'] = [this_antenna.nn_temp_var_zs[jd] for jd in jds]
z_score_cols['ee Temporal Discontinuties Modified Z-Score'] = [this_antenna.ee_temp_discon_zs[jd] for jd in jds]
z_score_cols['nn Temporal Discontinuties Modified Z-Score'] = [this_antenna.nn_temp_discon_zs[jd] for jd in jds]
for col in z_score_cols:
    df[col] = z_score_cols[col]

ant_metrics_cols = {}
ant_metrics_cols['Average Dead Ant Metric (Jee)'] = [this_antenna.Jee_dead_metrics[jd] for jd in jds]
ant_metrics_cols['Average Dead Ant Metric (Jnn)'] = [this_antenna.Jnn_dead_metrics[jd] for jd in jds]
ant_metrics_cols['Average Crossed Ant Metric'] = [this_antenna.crossed_metrics[jd] for jd in jds]
for col in ant_metrics_cols:
    df[col] = ant_metrics_cols[col]

redcal_cols = {}
redcal_cols['Median chi^2 Per Antenna (Jee)'] = [this_antenna.Jee_chisqs[jd] for jd in jds]
redcal_cols['Median chi^2 Per Antenna (Jnn)'] = [this_antenna.Jnn_chisqs[jd] for jd in jds]   
for col in redcal_cols:
    df[col] = redcal_cols[col]

# style dataframe
table = df.style.hide_index()\
          .applymap(lambda val: f'background-color: {status_colors[val]}' if val in status_colors else '', subset=['A Priori Status']) \
          .background_gradient(cmap='viridis', vmax=mean_round_modz_cut * 3, vmin=0, axis=None, subset=list(z_score_cols.keys())) \
          .background_gradient(cmap='bwr_r', vmin=dead_cut-.25, vmax=dead_cut+.25, axis=0, subset=list([col for col in ant_metrics_cols if 'dead' in col.lower()])) \
          .background_gradient(cmap='bwr_r', vmin=crossed_cut-.25, vmax=crossed_cut+.25, axis=0, subset=list([col for col in ant_metrics_cols if 'crossed' in col.lower()])) \
          .background_gradient(cmap='plasma', vmax=4, vmin=1, axis=None, subset=list(redcal_cols.keys())) \
          .applymap(lambda val: 'font-weight: bold' if val < dead_cut else '', subset=list([col for col in ant_metrics_cols if 'dead' in col.lower()])) \
          .applymap(lambda val: 'font-weight: bold' if val < crossed_cut else '', subset=list([col for col in ant_metrics_cols if 'crossed' in col.lower()])) \
          .applymap(lambda val: 'font-weight: bold' if val > mean_round_modz_cut else '', subset=list(z_score_cols.keys())) \
          .applymap(lambda val: 'color: red' if val > mean_round_modz_cut else '', subset=list(z_score_cols.keys())) \
          .bar(subset=list(bar_cols.keys()), vmin=0, vmax=1) \
          .format({col: '{:,.4f}'.format for col in z_score_cols}) \
          .format({col: '{:,.4f}'.format for col in ant_metrics_cols}) \
          .format('{:,.2%}', na_rep='-', subset=list(bar_cols.keys())) \
          .set_table_styles([dict(selector="th",props=[('max-width', f'70pt')])]) 

Table 1: Per-Night RTP Summary Info For This Atenna

This table reproduces each night's row for this antenna from the RTP Summary notebooks. For more info on the columns, see those notebooks, linked in the JD column.

display(HTML(f'<h2>Antenna {antenna}, Node {this_antenna.node}:</h2>'))
HTML(table.render(render_links=True, escape=False))

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
/tmp/ipykernel_4136/880101968.py in <module>
----> 1 display(HTML(f'<h2>Antenna {antenna}, Node {this_antenna.node}:</h2>'))
      2 HTML(table.render(render_links=True, escape=False))

AttributeError: 'NoneType' object has no attribute 'node'

Load antenna metric spectra and waterfalls from auto_metrics notebooks.

htmls_to_display = []
for am_html in auto_metric_htmls:
    html_to_display = ''
    # read html into a list of lines
    with open(am_html) as f:
        lines = f.readlines()
    
    # find section with this antenna's metric plots and add to html_to_display
    jd = [int(s) for s in re.split('_|\.', am_html) if s.isdigit()][-1]
    try:
        section_start_line = lines.index(f'<h2>Antenna {antenna}: {jd}</h2>\n')
    except ValueError:
        continue
    html_to_display += lines[section_start_line].replace(str(jd), f'<a href="{jd_to_auto_metrics_url(jd)}" target="_blank">{jd}</a>')
    for line in lines[section_start_line + 1:]:
        html_to_display += line
        if '<hr' in line:
            htmls_to_display.append(html_to_display)
            break

Figure 1: Antenna autocorrelation metric spectra and waterfalls.

These figures are reproduced from auto_metrics notebooks. For more info on the specific plots and metrics, see those notebooks (linked at the JD). The most recent 100 days (at most) are shown.

for i, html_to_display in enumerate(htmls_to_display):
    if i == 100:
        break
    display(HTML(html_to_display))

Antenna 137: 2459787

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	ee Shape Modified Z-Score	nn Shape Modified Z-Score	ee Power Modified Z-Score	nn Power Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Discontinuties Modified Z-Score	nn Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	48.282939	26.353169	28.385360	48.282272	48.282939	33.871536	33.444930	8.771705	10.619451

---

Antenna 137: 2459786

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	59.538503	27.724872	27.643152	59.538503	58.437418	50.088551	52.392285	16.520611	12.907417

---

Antenna 137: 2459785

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	59.346180	18.973044	18.015679	59.346180	57.829473	43.485570	45.224777	1.005280	0.700386

---

Antenna 137: 2459784

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	ee Shape Modified Z-Score	nn Shape Modified Z-Score	ee Power Modified Z-Score	nn Power Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Discontinuties Modified Z-Score	nn Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	59.544209	25.748767	26.068939	59.020282	59.544209	49.767751	47.406759	0.786363	0.720344

---

Antenna 137: 2459783

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Shape	24.758169	24.758169	23.386129	17.785466	16.942298	22.371822	24.347660	5.481641	4.855502

---

Antenna 137: 2459782

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Shape	22.500399	22.500399	21.617420	8.920546	8.328477	12.605915	13.737945	2.823113	2.711700

---

Antenna 137: 2459781

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	ee Shape Modified Z-Score	nn Shape Modified Z-Score	ee Power Modified Z-Score	nn Power Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Discontinuties Modified Z-Score	nn Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Temporal Variability	16.162880	3.262453	3.547004	6.437916	6.625634	16.118140	16.162880	1.185535	1.264598

---

Antenna 137: 2459778

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	ee Temporal Variability	23.749857	15.311558	14.337225	10.939860	10.283551	22.088167	23.749857	2.164046	1.796926

---

Antenna 137: 2459776

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	ee Shape Modified Z-Score	nn Shape Modified Z-Score	ee Power Modified Z-Score	nn Power Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Discontinuties Modified Z-Score	nn Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Temporal Variability	37.573030	21.343091	23.894450	7.913494	7.897085	36.298492	37.573030	1.241990	1.484427

---

Antenna 137: 2459773

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	80.361804	25.337733	21.305493	80.361804	76.857595	19.610479	20.024814	10.314406	7.801303

---

Antenna 137: 2459772

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	ee Shape Modified Z-Score	nn Shape Modified Z-Score	ee Power Modified Z-Score	nn Power Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Discontinuties Modified Z-Score	nn Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	ee Shape	nan	nan	nan	inf	inf	nan	nan	nan	nan

---

Antenna 137: 2459771

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	30.706836	15.142841	11.253106	30.706836	27.739521	3.911758	3.288611	4.798867	3.248270

---

Antenna 137: 2459770

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	52.930028	24.372233	8.585994	52.930028	34.956516	13.008543	10.981400	5.706936	3.955434

---

Antenna 137: 2459769

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	67.321593	27.335882	12.174089	67.321593	53.849485	16.645103	22.433915	5.242840	6.931968

---

Antenna 137: 2459768

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	ee Shape Modified Z-Score	nn Shape Modified Z-Score	ee Power Modified Z-Score	nn Power Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Discontinuties Modified Z-Score	nn Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	48.750419	15.016710	32.282850	33.577574	48.750419	13.925663	15.836910	10.785044	11.162352

---

Antenna 137: 2459767

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	ee Shape Modified Z-Score	nn Shape Modified Z-Score	ee Power Modified Z-Score	nn Power Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Discontinuties Modified Z-Score	nn Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	ee Shape	nan	nan	nan	inf	inf	nan	nan	nan	nan

---

Antenna 137: 2459766

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	30.391898	18.004522	6.787890	30.391898	16.389633	4.493630	18.402994	4.553127	6.367703

---

Antenna 137: 2459765

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	nn Shape Modified Z-Score	ee Shape Modified Z-Score	nn Power Modified Z-Score	ee Power Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Discontinuties Modified Z-Score	ee Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Shape	nan	nan	nan	inf	inf	nan	nan	nan	nan

---

Antenna 137: 2459763

Ant	Node	A Priori Status	Worst Metric	Worst Modified Z-Score	ee Shape Modified Z-Score	nn Shape Modified Z-Score	ee Power Modified Z-Score	nn Power Modified Z-Score	ee Temporal Variability Modified Z-Score	nn Temporal Variability Modified Z-Score	ee Temporal Discontinuties Modified Z-Score	nn Temporal Discontinuties Modified Z-Score
137	N07	RF_maintenance	nn Power	35.905878	9.618597	21.889488	4.915890	35.905878	26.874632	4.816486	17.813942	4.274940

---